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Abstract 

This paper presents and implements an iterative feedback design algorithm for stabilisation of 
discrete-time switched systems under arbitrary switching regimes. The algorithm seeks state feed- 
back gains so that the closed-loop switching system admits a common quadratic Lyapunov function 
(CQLF) and hence is uniformly globally exponentially stable. Although the feedback design problem 
considered can be solved directly via linear matrix inequalities (LMIs), direct application of LMIs for 
feedback design does not provide information on closed-loop system structure. In contrast, the feed- 
back matrices computed by the proposed algorithm assign closed-loop structure approximating that 
required to satisfy Lie-algebraic conditions that guarantee existence of a CQLF. The main contribu- 
tion of the paper is to provide, for single-input systems, a numerical implementation of the algorithm 
based on iterative approximate common eigenvector assignment, and to establish cases where such 
algorithm is guaranteed to succeed. We include pseudocode and a few numerical examples to illus- 
trate advantages and limitations of the proposed technique. 

1 Introduction 

We consider the discrete-time switching system (DTSS) 

= A m x k + B i(k) u k , (1) 

with x k e R n , u k G K. m , defined by a switching function 

i(k) G n := {1,2, ...,n}, for all k, 

and a set of controllable subsystem pairs {(A, B,) : ' G n}, where the input matrices B,-, /' g n, have full 
column rank. 

We address feedback control design of the form 

u k = K i{ k)X k , (2) 
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(which assumes that at every time instant k the "active" subsystem given by i(k) is known) so that the 
resulting closed-loop system 

x k+1 = A^jXfc, where (3) 
A ( CL = Aj + B,K,, for / e n, (4) 

be exponentially stable under arbitrary switching. 

It is well-known that ensuring that the closed-loop matrices A ( CL are stable for each / e n is necessary 
but not sufficient to ensure the stability of the DTSS ©-dD under arbitrary switching A necessary 
and sufficient condition for uniform exponential stability under arbitrary switching is the existence of 
a common Lyapunov function for each of the component subsystems in ©-(ID [2]. Such Lyapunov 
functions, however, will in general have complex level sets, which makes their numerical search difficult 

The search for common quadratic Lyapunov functions (CQLF), although restrictive, is appealing, 
since these functions play an important role in the stabilisation of linear time-invariant systems such as 
the component subsystems in Q}. The design of feedback matrices in © so that the DTSS ©-I© 
admits a CQLF may be posed as follows. 

Problem 1 . Given the matrices A, e R nxn and S, e R nxm for i e n, design feedback matrices K, e R mx " 
such that the DTSS closed-loop system 0)-® admits a CQLF. 

Quadratic Lyapunov functions are amenable to efficient numerical computation using linear matrix 
inequalities (LMIs). For example, Problem [Qcan be solved by finding X = X T > and W, to satisfy the 
LMIs 



X (AjX + BjNif 
A,X + BjNj X 



> 0, / € n, (5) 



where the required feedback matrices are given by K, = N,X~ ] and the CQLF is V(x) = x r X~ 1 x (see for 
example @|, 0]). An advantage of this approach is that the feasibility of the LMIs OD is necessary and 
sufficient for the DTSS considered to admit a CQLF. However, blind application of such control design 
strategy gives no insight on the structure of the closed-loop DTSS. Thus, as pointed out in Jg|, these LMI 
methods "lack transparency and interpretability that was a feature of classical techniques" and hence "a 
pressing need remains for analytic tools to support the design of stable switched systems" (also see @|, 
|7|| for similar comments). 

As an alternative to the LMI approaches, the authors in @] propose a pole-placement technique for 
single-input single-output continuous-time switching systems. The strategy in |§| seeks to guarantee 
closed-loop uniform global exponential stability under arbitrary switching by designing controllers that 
achieve a closed-loop common eigenvector structure. By constraining such eigenstructure and the class 
of controllers allowed, the strategy in |6|| simplifies the design process, providing analytically transparent 
solutions in a restricted but practically important class of systems. 

The present paper presents another strategy that seeks to "activate" analytic tools into a feedback 
design methodology (much in the spirit of [SQ) to solve Problem [Q Our strategy follows from the previous 
paper [SO, which introduced an iterative algorithm to seek feedback gains that make the set of closed- 
loop subsystems ©-© satisfy Lie algebraic conditions that guarantee the existence of a CQLF (To| . 
While such Lie-algebraic conditions are restrictive, since they are not necessary for the existence of a 
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CQLF, we believe they offer an insightful way to understand and exploit fundamental system structure in 
feedback design for DTSS. 

The theoretical results in [9] brought forward the following important consequences: (i) if the pro- 
posed Lie-algebraic feedback design problem is feasible, its solution can be found in an iterative manner 
(similarly to the way solvability of Lie-algebras can be checked for autonomous switched systems) by 
seeking feedback gains that assign a single common stable eigenvector at each iteration step; and (ii) 
in seeking such common eigenvector, if at any iteration step more than one vector can be assigned by 
feedback, it is irrelevant which one is chosen by the procedure. These observations provide motivation 
to the development of numerical implementations, which were not discussed in @|. 

The present paper focuses on the numerical implementation aspects of the iterative algorithm in- 
troduced in |9||. A key question in the proposed approach is the lack of robustness of the Lie-algebraic 
conditions that are sought to satisfy by feedback design. Indeed, it is well-known that these conditions 
are destroyed by arbitrarily small perturbations to the individual matrices [1o[ §2.2.4], and thus, there are 
a priori no guarantees that the algorithm in [9J] can find any solution at all in a (necessarily approximate) 
numerical implementation. 

However, suppose that for a given set of systems to be stabilised there exists, in a neighbourhood of 
the original, a feasible (a priori unknown) set of systems, which is stabilisable and such that the resulting 
closed-loop systems satisfy the Lie-algebraic conditions that guarantee the existence of a CQLF. Then, 
by continuity of such CQLF, and if the neighbourhood is sufficiently small, there will also exist a CQLF 
for the original set of systems, despite the fact that it may be impossible to make them satisfy the 
Lie-algebraic conditions. 

The main contribution of the present paper is to, based on the above argument, derive and math- 
ematically justify a specific numerical implementation, for single-input systems, of the design algorithm 
proposed in §| that will succeed not only in cases where the Lie-algebraic conditions are satisfied but 
also in approximate cases. The proposed numerical implementation is based on the solution of an 
optimisation problem that seeks feedback matrices that will achieve closed-loop systems with an ap- 
proximate common eigenvector. The existence of a CQLF for the corresponding closed-loop DTSS may 
be readily checked a posteriori with a set of informed LMIs. This step is necessary, since it is in general 
a priori unknown if the "exact" problem is feasible in some neighbourhood of the original system data. 
The algorithm has been implemented in Matlab®, and numerical examples to illustrate its application 
and discuss its advantages and limitations are provided. 

The rest of the paper is organised as follows. The proposed algorithm and its core, a procedure for 
approximate common eigenvector assignment (Procedure CEA), are presented and explained in Sec- 
tion|2] The main theoretical results, justifying the numerical implementation of the proposed algorithm for 
single-input systems, appear in Section [3] Section g] presents details about the Matlab® implementation 
of the algorithm and some illustrative numerical examples, and Section [5] gives the paper conclusions. 
The main proofs (Theorems [2] and [3} are given in the Appendix. 

Notation R and C denote the real and complex numbers. || • || denotes Euclidean norm or correspond- 
ing induced matrix norm. M* denotes the conjugate transpose of M. p(-) denotes spectral radius and 
(M); ik is the k-th column of M. If M e C" xm , Im M denotes {x e C" : x = My,y e C m }. The Euclidean 
distance between a vector i/eC" and a set V c C" is denoted d(v, V). j denotes v/^T. 
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2 Feedback Control Design 

A sufficient condition for the closed-loop DTSS |[3)-||4]| to admit a CQLF is given by the following result, 
which is a minor modification of [Til Theorem 6.18]. 

Lemma 1 . If p(Af L ) < 1 for i e n, and the Lie algebra generated by {Af L : i e n} is solvable, then ©-© 
admits a CQLF. 

In matrix terms, the fact that the Lie algebra generated by {,4 ( CL : / e n} be solvable is equivalent to 
the existence of a single invertible matrix T e c nx " such that T~M ( CL T is upper triangular for / e n (even 
if /A ( CL have real entries, those of T may be complex [12]). 

2.1 The algorithm 

In 0], we established that given A, and B,, there exist feedback matrices K, that cause the Lie algebra 
generated by /\ ( CL to be solvable if and only if such feedback matrices can be computed by an algorithm 
that performs iterative common eigenvector assignment by feedback. A matrix version of such algorithm 
is given in pseudocode as Algorithm [TJ 

At every iteration, Algorithm [TJ executes Procedure CEA [see ©]. This procedure attempts to find 
feedback matrices Ff and a vector v\ so that (Af + BfFf)v^ = A^vf and |Af| < 1 for / e n, i.e. so 
that v\ becomes a common eigenvector of a set of closed-loop matrices, with corresponding stable 
eigenvalues. The parameters e c and given as arguments to Procedure CEA are required for numerical 
reasons, and will be explained in Section [2T21 
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Algorithm 1 : Approximate triangularisation by feedback 



Data: A, e R nxn , B, e R nx1 for / e n, e c > 0, e d > 
Output: (J, K, for / e n 
begin Initialisation 

A; := Aj, Bf := B,, K° := 0, LA| := I ; 
U :=[] (empty), £^0; 

repeat 

£ <- t + 1 , n r ^ n - £ + 1 
tf. {Ff <" CEA({Af , B/}^ , e C! e d ) (6) 



Define 

A/< CL := A/ + Bf Ff, (7) 

r=£ 

Kf <- K/- 1 + Ff(f] U;) (8) 

r=1 

if £ < n then 

Construct unitary matrix with first column v\: 

[vi\4\---\vf lr ]eC n ' xn '. (9) 
Assign 

u M ^[4\ ■■■Ki 00) 

A/ +1 ^U| +1 Af CL U, +1] (11) 
Bf + '^U* M Bf, (12) 

until £ = n; 



Remark 1 . /f is straightforward to check that if Algorithm\1}terminates successfully and at every iteration 
(£ = 1 , . . . , n) Procedure CEA is able to find the vector v\ and feedback matrices Ff such that (Af + 
BfFf)v\ = for i £ n, then the matrices A ( CL given by © with K, as computed by the algorithm are 
such that U*Af L U are upper triangular and is the i-th main-diagonal entry of U*Af L U. 

Note that a slight modification of Algorithm [Q is necessary to ensure that real feedback matrices K, 
are computed. We do not explain such modification here due to space limitations, and because it does 
not add any essential information to our main results. The implemented computational version of the 
algorithm, employed in Section m does indeed ensure such condition. 

Our result in [9] established that a vector v\ and feedback matrices Ff such that {Af + BfFf)v\ = \ l .v\ 
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with | | < 1 , for / G n, will exist at every iteration of Algorithm Q] if and only if there exist K, such that A ( CL 
as in d3) generate a solvable Lie algebra and satisfy p(A ( CL ) < 1 . Such result was of a theoretical nature, 
since in a numerical implementation determining whether a vector is exactly an eigenvector or is close 
to being so is an extremely difficult task |13|. 



2.2 Approximate common eigenvector assignment 

In this section, we provide a numerical implementation of Procedure CEA executed by Algorithm [Hand 
establish some of its properties. We give Procedure CEA in pseudocode first, and next define and 
explain each of its parts. 



Procedure CEA 

Input: Af e C" rX \ Bf e C" rX1 , for / e n, e c , e d 
Output: v\, Ff for / e n 
if n r = 1 then 

Select Af< CL such that |A £ ' CL | < 1 - e c ; 

vf ^ 1 , Ff <- -(ef)-Mf + A^ ■ 

else 

if S(e c ,e d ) = 0then 
| Stop: unsuccessful termination, 
else 

vf <- argmin v&S(eoied) J(v) ; 
_ Ff ^ M,(vf) ; 



Note that n r is the state dimension, which decreases by 1 at every iteration of Algorithm [Q Hence, 
the case n r = 1 in Procedure CEA corresponds to the trivial case of a one-dimensional single-input 
system. If every subsystem is controllable, then + Bf e C and hence Af' CL can be arbitrarily chosen. 

For the case n r > 1 , Procedure CEA utilises the matrices 



EM 
H,(v) 
M,(v) 
Af< c \v) 



= (vv* - \)Af, 
= (w* - \)Bf, 

= -(H(v)*H / (v))- 1 H / (v)*E / (v), 
= Af + BfMi(v), 



(13) 
(14) 
(15) 
(16) 



the sets 



5i 
S 2 (e c ) 



= {v e C" r : ||v|| = 1} 



■<n, 



\Af' CL (v)v\\ < 1 - e £ 



N 



5(e c ,e d ) := S^ n 5 2 (e c ) n 5 3 (e d ), 



(17) 
(18) 

(19) 
(20) 
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and the cost function 



J(v) := J \\[Ei(v) + H,(v)M,(v)]v 



(21) 



If n r > 1 , Procedure CEA checks whether the set S(e c , e d ) given by d^UJ is empty. Comments on the 
case <S(e c , e^) = are given later in Remark[2] If S(e c , e d ) + 0, Procedure CEA searches for a vector v\ 
that minimises J(v) as given by flU), subject to the constraint v e <S(e c , e d ). The constraint set <S(e c , e d ) 
in (|2"0t is the intersection of three sets: 5i in fT7). which constrains the search to unit vectors; S2(e c ) in 
(QD, which imposes the stability constraint |Af | < 1 , as we will shortly demonstrate; and S 3 (e d ) in ff9) . 
which is included for technical reasons discussed next and justified in the next section. 

Note that if e d > and v e S 3 (e d ), then (T9]| implies that v £ Im , for / e n. It follows that /-/,■( v) 
has the same (column) rank as Bf [see (fT4)] and, in this case, Mj(v) in COD is well-defined if Bf has full 
column rank (l.i. columns). 

Thus, Procedure CEA requires Bf to have full column rank. At the first iteration of Algorithm [Q 
we have B- = B, and hence B- has full column rank by assumption. At subsequent iterations of the 
algorithm, such condition is ensured by the following result. 

Lemma 2. Let Bf have full column rank and let v\ be computed by Procedure CEA. Then, Bf +] has full 
column rank. 

Proof. Since v\ is effectively computed by Procedure CEA and since S(e c , e d ) is topological^ closed, 
then v\ G <S(e c , e d ) and hence v\ £ Im Bf. The result then follows straightforwardly from the fact that the 
columns of © form a basis, (QO) and C[2). □ 

We are now ready to show in what sense the minimisation of J(v) is related to the assignment of a 
common eigenvector. 

Lemma 3. Let e c > 0, e d > and v e <S(e c , e d ). Then, 
i) J(v) > 0, 

//; J(v) = if and only if Af' CL (v)v = Afv with |Af| < 1 - e c . 
///; There exists G, e C mx " r such fhaf 



(A/ + Bf6,)v = Afv, 



(22) 



/fand on/y/f/\f' 0L (v)v = Afv. 

Proof, [i) Straightforward from {21} and since J is well defined on <S(e c , e d ). 
(=>) J(v) = implies that 



= [EM + H,{v)M,{v)]v 
= \\(vv* - \)[Af + BfMi(v)]v 



(23) 



for / G n, where we have used (13), (JT4) and ||2T}. In turn, ||23| implies that 



[Af + BfMi(v)]v = Af- CL (v)v = Xjv, 



(24) 
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where we have used (TB). From (2J), |Af | = ||A (v) v|| < 1 — e c < 1 since v e <S(e c , e d ). 
(<=) Left-multiply (24) by (vv* - I) to obtain (25), which implies J(v) = 0. 

ED (=>) Left-multiplying (22) by (vv* - I) yields [£,(v) + H ; (v)G,]v = 0. Let u-, := G,v. Then, E,(v)v = 
— Hj(v)Uj. Since H,(v) has full column rank, then u-, = Mj(v)v = G,v and hence by (22) A e,CL (v)v = 
[Af + Bf Mi(v)]v = Af v. 

(<H Just take G; = M,(v). □ 

Lemma[5J shows that, for v e S(e c ,e d ), J(v) = if and only if v can be assigned by feedback as 
a common eigenvector with corresponding stable eigenvalues, and J(v) > otherwise. Therefore, the 
search for a vector that minimises J(v) can be interpreted as the search for a vector that is closest to an 
assignable common eigenvector. 

Remark 2. If S(e c , e d ) = 0, then no feasible vector exists for the optimisation performed by Proce- 
dure CEA. Therefore, no vector can be computed, and Procedure CEA and hence Algorithm U} will 
terminate unsuccessfully. In the next section, we show that there exist e c , e d > so that S(e c , e d ) + in 
the cases of interest. 

In the next section, we establish conditions under which the vector and feedback matrices computed 
by Procedure CEA cause Algorithm [Qto yield feedback matrices so that the closed-loop DTSS admits 
a CQLF. 



3 Main Results 

We now derive several results that justify the use of Algorithm [T] for feedback control design of single- 
input systems. Recall that the constraint set <S(e c , ed) in the optimisation solved in Procedure CEA 
forces the search to be performed over vectors v that satisfy v <£ Im Bf. The following result justifies this 
constraint for single-input DTSS with controllable subsystems. 

Lemma 4. Lef n r > 1 , Af € £ n,xn ' , Bf eC"' x1 , (Af, Bf) be controllable and suppose that v + Qand Ff 
satisfy (Af + BfFf)v = A,v, for i £ n. Then v <£ Im Bf for i e n. 

Proof. Suppose for a contradiction that v £ Im Bf for some k e n. Then v = Bfu for some + u £ C. 
We have (Af + BfFf)Bfu = XkBfu and hence Bf and AfBf are linearly dependent, which prevents the 
pair (Af, Bf) from being controllable because n r > 1 . □ 

Lemma |U establishes that, for a single-input DTSS with controllable subsystems, every vector that 
can be made a common eigenvector by feedback will not be contained in the image of any input matrix 
Bf. Therefore, in this case such constraint on the optimisation problem is justified. 

In the sequel, we will say that a set of matrix pairs is controllable if every matrix pair in the set is 
controllable. We also require the following definitions. 

Definition 1 (CEAS, 7-CEAS). A set of matrix pairs Z e = {(Af e C" rX " r , Bf e C" rX1 ) : / e n} is said 
to be CEAS (Common Eigenvector Assignable with Stability) if there exist + v e C" r , Af e C with 
| Af | < 1 , and matrices Ff e C 1 x " r , such that 

(Af + BfFf)v = Af v, for i e n, (25) 
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If Z e is CEAS, we say that v e and F-, e C 1x " r are compatible with Z l whenever (25$ holds for 
| Af | < 1 . If (25) holds with | Xj\ < 1 - 7 for some < 7 < 1 , we say that Z l is j-CEAS and refer to the 
corresponding v and Ff as 7 -compatible with Z e . 

Definition 2 (SLASF, 7-SLASF). A set of matrix pairs Z = {(A; G C nxn , B, e C nx1 ) : / G n} /s sa/d to 
be SLASF (Solvable Lie Algebra with Stability by Feedback) if there exist K, G C 1x " such that Af 1 as 
/n 0) generate a solvable Lie algebra and satisfy p(Af L ) < 1 . If Z is SLASF, we say that Kj e C 1xn 
are compatible with Z if Af L as in 0) generate a solvable Lie algebra and satisfy p(A ( CL ) < 1 . If K\ exist 
so that, in addition, p(A ( CL ) < 1 - 7 for some < 7 < 1 , we say that Z is 7-SLASF and refer to the 
corresponding K, as 7 -compatible with Z. 

We next state a version of our previous result of @| for the specific case of single-input systems as 
follows. 

Theorem 1. Let Z = {(A, e C nxn ,B ( - e C nx1 ) : / e n}, consider AlgorithmWor some suitable choice 
of e e c and e e d , and let Z l = {{Af, Bf) : ;' g n} and < 7 < 1 . Then, Z /s 7-SLASF /f and on/y /f /s 
7-CEAS and v(, Ff are ^-compatible with Z l for l=\,...,n. 

Suppose that Z 1 = {(Af, Bf) : i e n} is CEAS. Whenever Ff, v\ are compatible with Z l and a unitary 
matrix having v\ as its first column is given: 

M\4\---\K- M l (26) 

then Af' CL , M , Af +: and Bf +1 will denote the matrices given by 0, (TO), 03} and (T2), respectively, 
when the hatted matrices are employed. 

The main result of this section, namely Theorem [3] requires the following preliminary theorem. 

Theorem 2. Lef Z l = {(Af, Bf) : / e n} be 7-CEAS and controllable, with Af e C n ' xn ' , Bf e C" rX1 , 
n r > 1 , < 7 < 1 . Consider the following sets 

7f a = {v G 5i : 3F e C 1 x " r , A e C such fhaf (Af + Bf F)v = Xv with |A| < 1 - 7}, (27) 

N 

Xy £ =n^5r' (28) 
/=1 

and fhe following quantity 

e e f := inf^m/n )eM d(v,lmBf). (29) 

Then, e 1 ^ > and each < 4 < 7 a" ' < e d < e d* ensure that for every e > fhere ex/'sfs a 
corresponding 5 > so fhaf for each Af , Bf satisfying 

\Af-Af\\<5, 
\Bf - Bf || < 5, 

fhere ex/sf vf e •S(e e c /2, e e d /2) and Ff compatible with Z l (\r\, Ff may depend on the specific Af, Bf), 
and a unitary matrix (25$ that cause 
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n -J fn 

|| v\ - v||| < e, 

\\F l - F e \\ < e 

where and Ff are the output of Procedure CEA with Af, Bf, e c = e e c , and = e e d as inputs. 



(30) 
(31) 




(32) 
(33) 



where Ue+] , Af^ and Bf +: are the matrices computed at iteration £ of Algorithm[1}from Af and Bf, 
with v\ and Ff as above. 

Theorem [2][i) establishes that if the matrices Af and Bf given as inputs to Procedure CEA are 
sufficiently close to some Af and Bf which form a CEAS set Z e , then the vector v\ and feedback 
matrices Ff computed by such procedure will be as close as desired to some v\ and Ff compatible 
with Z L . In general, whether the (given) matrices Af, Bf are sufficiently close to some Af, Bf with the 
required property will not be known. However, the significance of this result lies precisely in the fact that 
it establishes a type of continuity relation between the result of the procedure and an "exact" result, even 
if the latter result is not known. In addition, such continuity justifies the numerical implementation of the 
procedure, since numerical computation will always yield an approximate result. 

In broad terms, Theorem EM) shows that if, at step I of Algorithm [Q the matrices Af and Bf are 
sufficiently close to some "exact" ones, then the same will happen at step t + 1 . The proof of Theorem [2] 
is highly non-trivial and given in Appendix IA. 11 

We are now ready to state the main result of the paper. 

Theorem 3. Let Z = {(A; e C nx ", 6, e C nx1 ) : / e n} be SLASF and controllable. Then, there exist 
e*, e* > such that each < e c < e* and < e d < e* ensure that 

i) For every e > fhere exists a corresponding S > so that for each A,, B, satisfying 



there exist Kj compatible with Z (K, may depend on the specific A,, B,) that cause p(Af L ) < 1 - e c /2 
and 



where Af L = A; + B,K;, and Af 1 satisfies 0) with Kj obtained as output of Algorithm\1\with Aj, B h e c , 
and e d as inputs. 



A, - A|| < S, 
Bi - B/|| < 5, 



(34) 
(35) 



Af L 



< e 



(36) 



ii) There exists e > for which the closed-loop DTSS admits a CQLF, provided {34])-{35|) are 

satisfied with 6 corresponding to e as /n0) above. 
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Theorem [3] establishes that Algorithm [TJ will compute suitable feedback matrices not only in the 
"exact" case when the given A,, B, form a SLASF set, but also when they are close to other (possibly un- 
known) matrices A„ 8, with such property. This justifies the use of Algorithm [TJfor control design, since 
it gives a kind of robustness result for the feedback matrices computed by the algorithm. Theorem [3] 
also establishes that suitable feedback matrices will be computed for all positive e c and respectively 
less than e* and e* d . The latter quantities may be not known, since they depend on the possibly unknown 
A, and B,-. Consequently, in theory the parameters e c and e d should be selected as small as computa- 
tionally possible. In practice, however, there is a tradeoff in the selection of e c since the smaller e c , the 
higher the chances of Algorithm [Qyielding unsuitable feedback matrices when more than one common 
eigenvector could be assigned by feedback with one of them having unstable corresponding eigenval- 
ues and another having stable ones (Procedure CEA could select a local optimiser at the boundary of 
the constraint set instead of a global optimiser in its interior). 

After executing Algorithm [TJto compute feedback matrices K,-, we can check whether the closed-loop 
DTSS ©-(|4) admits a CQLF by solving the following LMIs 

P = P T > 0, P - (Af L ) T PAf L > 0, for / G n. (37) 

Note that, in this case, LMIs are used only to check whether a CQLF for the closed-loop system exists, 
and not for feedback design. If these LMIs are feasible, then not only the closed-loop DTSS admits a 
CQLF but also we have structural information on the DTSS since K, are such that the /\ ( CL are suitably 
close to being simultaneously triangularisable. 



4 Examples 

We next provide some numerical examples to illustrate the advantages and limitations of the proposed 
feedback design strategy. For the numerical implementation of Procedure CEA, a feasible vector is first 
sought using Matlab® Optimization Toolbox function fgoalattain. If such vector is found, it is passed as 
initial point to the optimisation, implemented via the function fmincon. 



4.1 Randomly created DTSS 



The following subsystems were created randomly but such that p(A\) < 1 and p(A 2 1 ) < 1 . 



A, = 
A 2 = 



0.574 0.074 0.089 1 

0.074 0.572 -0.091 , Bi = 

0.089 -0.091 0.538 J 

-0.737 0.386 -1.680 
1.351 0.638 0.035 , B 2 = 
1.071 -1.295 -0.936 



-0.038 
0.327 
0.175 



0.114 
1.067 



(38) 
(39) 



Executing Algorithm [TJ choosing e c = e d = 10 4 yields a feasible optimisation at every iteration and 
returns 



Ki =[ -3.6480 -7.2304 8.7751 ], / / _ 
K 2 =[-0.3159 2.0235 0.2695], U ~ 



0.4647 0.8287 -0.3120 
-0.7770 0.2126 -0.5925 
-0.4246 0.5178 0.7427 



(40) 



It can be shown that LMIs (37) are feasible and hence the closed-loop DTSS ©-0 with K, as in (40) 
admits a CQLF and is hence stable under arbitrary switching. In this case, it can be checked that 
U*Af-U are not upper triangular but are close in the sense that the entries below the main diagonal are 
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small. Therefore, the use of feedback matrices (HU) designed via Algorithm [JJprovides some insight into 
the structure of the closed-loop DTSS. 

On the other hand, solution of the LMIs © for both feedback design and CQLF computation, which 
yields 

K-\ =[-1.2267 -0.7211 -1.8731 ], K% = [ -0.5140 1.3826 1.1613] , 

is guaranteed to produce a closed-loop DTSS stable under arbitrary switching but provides no structural 
information. 



4.2 DTSS with no CQLF 

Consider the systems 

= [?]• 



For a = 1.5, the LMIs © are not feasible. Therefore, no CQLF exists for this DTSS. Executing Algo- 
rithm [TJ with e c = e d = 10 -4 yields an infeasible optimisation at the first iteration, S{e c , e^) = 0. 

For a = 1 .4999, the LMIs © are feasible. In this case, Algorithm [TJfor the selected e c , ed yields an 
infeasible optimisation. However, reducing e c to 10~ 5 allows the algorithm again to compute suitable 
feedback matrices. 



4.3 CQLF exists but Algorithm Q] fails 

Consider again the DTSS <38t-|[39l, with the addition of the subsystem 
As = 



0.352 0.159 


-1.129" 


, B 3 = 


' -0.433 


0.159 


0.262 





-1.129 0.262 


-0.705. 








Algorithm [TJfor e c = ed = 10 4 yields 

Ki =[-15.3542 3.8969 -11.3814], 
K 2 = [ 0.0734 0.9747 2.7288 ] , 
K 3 = [ -1.3542 0.8334 -4.5001 ] , 



(41) 



U = 



0.1662+0.0234/ -0.0918-0.6508/ -0.7347 
0.9744+0.1374/ 0.0093+0.0662/ 0.1650 
0.0587+0.0083/ 0.1049+0.7433/ -0.6580 



However, the optimisation informs that there is an active inequality, corresponding to the stability con- 
straint |T8). In this case, the LMIs (|37) are not feasible and hence no CQLF exists when the feedback 
matrices dJTJ are employed. 

On the other hand, the LMIs © are feasible and hence other feedback matrices may indeed produce 
a closed-loop DTSS with a CQLF. 



5 Conclusions 

This paper complements the theoretical results in @, and contributes to furthering the understanding 
of fundamental system structure in feedback stabilisation of DTSS. We have presented a numerical 
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implementation of a feedback design strategy for DTSS based on Lie-algebraic solvability. The proposed 
strategy seeks feedback matrices to achieve a closed-loop system structure that approximates that 
required to satisfy such Lie-algebraic stability criteria. 

The main theoretical contribution of the paper establishes that if a system for which the Lie-algebraic 
conditions considered exists in a suitably small neighbourhood of the given system data, then our im- 
plementation will find feedback matrices so that the corresponding closed-loop DTSS admits a CQLF 
even if the considered Lie-algebraic conditions are not met by the given system. However, since the 
existence of such feasible "exact" system is in general unknown, the resulting closed-loop system is 
not guaranteed to admit a CQLF but the latter may be checked with a set of informed LMIs built with 
the computed closed-loop matrices. Whether such "exact" system exists suitably close to a given sys- 
tem has not been discussed, and remains a topic for further research. Future work will also consider 
extensions to multiple input systems, which are nontrivial and will possibly require the consideration of 
controllability indices in the algorithm. 



A Appendix 

A.1 Proof of Theorem |2] 

Throughout this proof, a hatted expression denotes the expression given by the corresponding equations 
when matrices Af, Bf are substituted for Af, Bf. For example, fij(v) = (vv* - \)Bf. 

We begin by establishing that e 1 ^ > 0. The set Ty is the set of all unit vectors that are 7-compatible 
with Z e , and since Z l is 7-CEAS, then Ty + 0. Since Z 1 is controllable, then (Af, Bf) is controllable and 
hence Bf + 0. 

Claim 1 . The set TJf is compact. 

Proof. The set ff is bounded since ff c <Si by definition. The set Tf can be equivalents defined as 
jf n = {veSi :3XeC such that Pf(X\ - Af)v = 0, with |A| < 1 - 7}, (42) 
where we have defined 



\-Bf((&f)*Bfy' (Bfy 



(43) 

Note that since + Bf e C" rX1 , then Pf is well-defined. Consider a sequence {v k }^ Q such that 
v k 6 7~ £ 7 for all k > and Hindoo v k = v. Since || Vfe [| = 1 for all k > and by the continuity of norms, 
then ||v|| = 1 necessarily. For every k > 0, we have 

Pf(X k \-Af)v k = 0, (44) 

for some A^eC with \X k \ < 1 — 7. From ((44) . 

lim Pfv k X k = lim \pf(v k - v)X k + PfvX k ] = PfAfv. (45) 

Since X k is bounded and lim/^oo v k = v, it follows that 

lim PfvX k = PfAfv. (46) 

k— >oo 
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If Pfv t 0, then Wm^^ X k = A with |A| < 1 - 7. If Pfv = 0, then PfAfv = 0. In either case, v G if 
and hence if is closed. Therefore, Zjf is closed since it is the intersection of a finite number of closed 
sets. □ 

The quantity e d * as defined in {29} is the infimum, over all vectors 7-compatible with Z e , of the 
minimum of the distance between such vectors and ImBf. Since (Af, Bf) is controllable and n r > 1, 
then Lemma 0] implies that v <f Im Bf for / e n and every v G 7^f. Therefore, for every v e TI, we have 
min/Qj d(v, Im Bf) > because Im Bf is closed for / G n. Since T~ is compact and min, eN d(v, Im Bf) 
is continuous and positive at every v G J~Jj, it follows that min/e, d(v, Im Bf) achieves a minimum on Zjf 
and hence e^'* > 0. 



A.1.1 Proof of Theorem HID 

Claim 2. For every < 4 < 7 < e e d < e £ f, there exists 5 2 > such friaf if n <S(e£, 4) * for a// 
Af , Bf satisfying 

\\Af-Af\\ < 62, \\Bf - Bf\\ < 5 2 . (47) 

Proof. Let and G, be 7-compatible with ^ and ||^|| = 1. Note that $ G Zf. By (US), G, and ^ 
satisfy, for / g n, 

(Af + BfG,y = A,^, (48) 

|A/|<1- 7 <1, (49) 

d(v\ Im Bf) > ef > 0. (50) 

Consider H-,(\fi) and /-//(v 11 ) from fi"4} . By ((50) and since Bf have full rank, then Hj(v^) has full rank. 
Then, for all 5o > sufficiently small, Hj(v^) has full rank whenever ||Bf - Bf\\ < 5q. Whenever the latter 
holds, the expression fi"6) is continuous on the entries of Af and Bf. Then, given e > 0, we can find 
< <5i < S such that 

\\Af' CL (v^ - Af- CL (v^\\ < e 
(WA 6 - AH < 81 

whenever jjj^ _ ^jj ^ < ^' (51) 

By |gE} and LemmalMn)- we have Af CL (v^)v^ = A/vl Since ||v^|| = 1, then 
||Af'°V)^ll = I A,-| < 1 - 7 < 1 

by d49}, and we can select e > small enough so that ||A ( . ,0L (v | i)v' , || < 1 — e e c and hence g 52(4) 
whenever (|5T} holds. For < a < 61, define 

d(a) := min inf d(v\\m Bf) 

i& Bf:\\&f-Bf\\<a 

and note that, by the continuity of d(v^, Im Bf) on the entries of Bf whenever Bf has full rank, and since 
by |50} d(0) > e £ f > 0, there exists 5 2 > sufficiently small for which 6(5 2 ) > e d > 0. Therefore, for 
such 5 2 we have $ G 5 3 (e^) and hence G S(e c , e d ) whenever (07} holds. □ 
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Claim 3. Consider < e e c < 7 and < e e d < . Then, there exists S 3 > so that 
i) El-El hold for every Af, Bf satisfying El). 

sup max \\Af' CL (v)v\\ < 1 - 4/2, (52) 
ves{44) ^ 

inf min d(v, Im Bf) > ei/2, (53) 

llA—A'll < s 3> \\Bf - Bf\\ < 5 3 . (54) 



//J For every e 2 > 0, friere exists < £4 < S 3 such that El ho/ds for every Af, Bf satisfying f56|) . 

|J(v) - J(v)| < e 2 , for all v G S(e £ c , e £ d ) (55) 
||Af-Af|| <<5 4 , l|Sf- ef|| < 5 4 . (56) 



Proof, [i) Consider the functions 

f({Af, Bf : i G n}) := sup max \\Af' 0L {v)v - Af' CL {v)v\\, 

g({Af , Bf : i g n}) := sup max \d(v, Im B/) - d(v, Im B/)|, 

which are non-negative whenever well-defined. Note that 

f({4f , Bf : / g n}) = = g({4f , Bf : / g n}), 

and that, since Bf full column rank, f, g are continuous on the entries of Af, Bf and well-defined when- 
ever d54t holds for some 63 sufficiently small. Therefore, for every e-| > 0, we can find a corresponding 
5 3 > so that 

f({< , Bf : / G n}) < „ , g({/\f , Bf : / G n}) < e, 

whenever (|54) holds. Select e-i < min{ec/2, e^/2}, and take the corresponding 5 3 . Thus, (I52"t-(l53"t hold 
whenever i(54) is true, 
[n) Consider the function 

h({Af,Bf:ieN}):= sup |J(v) - J(v)\. 

veS(ei,e l d ) 

Note that, because of (53), ri is well-defined, non-negative and continuous whenever <(54) holds and that 

h({Af, Bf : / G n}) = 0. 

Therefore, given e 2 > we can find < <5 4 < 5 3 , such that $55} holds whenever fSB) is true. □ 
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Fix e c and e e d so that < e e c < 7 and < e e d < . Consider Af and Bf satisfying both l@7) with 
5 2 as given by Claim |2] and l[5D with 63 from Claim ^G). Let v\ e S(e e c ,e d ) be a (global) minimiser of 
J(v) subject to v e S(e l c , e d ) [such minimiser exists because J is continuous on S{e c , e e d ) and S{e c , e d ) is 
compact and nonempty]. Define i c := e c /2 and e d := e^/2, pick e S(e c , e e d ) n Tjf, and note that 



< e c < 1 - maxmax(||A^ CL (vf)vfll. HA^O^Ill and 
< e d < min min {d(vf , Im Bf), d(v\ Im Bf)}, 



by <[52j>— <[53j) from Claim [M}, and since v\ G S(e c , e e d ) and G 5(4. ^ d )- Thus, we have vf G <S(e c , e d ), 
G <S(e c ,e d ) and by Lemma [MD. we know J(v^) = 0. Since v\ is a minimiser within S(e c ,e d ), and 
G <S(e^, e e d ), then J(vf) < J(v^). Thus, using Claim [3KID, we have 

J(vf) < J(&) = \3(&) - J(v*)\ < e 2 , (57) 
whenever \\Af - Af\\ < 5 A , \\Bf - Bf\\ < 64. (58) 

Also, using ((57) and Claim [3113, under ((58) we will have 

J(vf)<p(vf)-J(vf)|+J(vf)<2e 2 . 

Define V := G <S(e c , e d ) n <S(4, e d ) : J(v) = 0} and note that \fi G Vo- For each e 3 > 0, consider 

B(e 3 ) := {v G 5(e c , e d ) n 5(e c , e«) : d(v, V ) < e 3 }- 

By continuity of J and since both <S(e c , e d ) and 5(e c , e d ) are compact, for every e 3 > we can find e 2 > 
such that 

J(vf)<2e 2 =► v\eB{e z ) (59) 

Note that the closedness of <S(e c , e d ) and S(e e c , e l d ), and hence of S{e c , e d ) n <S(4> e^), is key in allowing 
the implication ((59). v\ G £>(e 3 ) implies the existence of v\ G Vo such that \\v\ - v\\\ < e 3 , and 
\t\ G Vo implies that J{v\) = 0. By Lemma[3]and since v\ G S(i c , id), we have Af ,cl {v\)v\ = a,v\ with 
|«/| < 1 — e c < 1 , which establishes that v\ and Pf := Mj(v^) are compatible with Z e . 

Next, since H,(v) has full rank whenever v G S(e c , e d )nS(ec, e^) and (58), then /W,(v) is continuous on 
v and on the entries of Af, Bf. Therefore, given e > we can find 5 5 > such that ||M/(vf) - M;(vf)|| < e 
whenever 

\\Af - Af\\ < S 5 , \\Bf - Bf\\ < 5 5 , || v>f - vf || < * 5 . 

Take e 3 = min{e, 5 5 } to select 64 as above. Finally, taking 5 = min{<5 2 , 64, 5 5 } concludes the proof of 
Theorem l2li). 

A.1 .2 Proof of Theorem |2"UI) 

Since (TTJ is continuous on the entries of U M and Af' CL , given e > we can find < 61 < e such that 
(|3"2l holds whenever 

\\Af* - Af*\\ < 6 U (60) 
\\0 M -U M \\<6,. 
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Since O is continuous on the entries of Bf and U M , then given e > we can find < S 2 < £1 so that 
holds whenever 

\\Bf-Bf\\ <5 2 , 

\\0 M - U M \\ < 6 2 . (61) 

Consider square unitary matrices W = [v^\U M ] and W = Note that, given v\, v\ and U M 

so that W is square and unitary, and ||vf || = 1 , for every 62 > we can find Um and < £3 < 62 so 
that 

\\W-W\\< 5 2 (62) 

whenever \\vf- - v^\\ < 5 3 . Note that |[62) implies ((61) since L/^+i and Um are the last n — £ columns of 
W and W, respectively. Since is continuous, given <5i > we can find < 64 < 5 3 so that $60$ holds 
whenever 



\Af - Af\\ < <5 4 , || Sf - 6/ 1| < 5 4 , and 

-£ 



\\Ff-Ff\\ < 6 4 . (63) 

Applying Theorem [HQ yields that each < e e c < 7, < e e d < e l f ensure that for the given 84 > 0, 
we can find < S < 64 so that for each Af, Bf satisfying \\Af - Af\\ < 5 and \\Bf - Bf\\ < 5, there exist 
vf G £(4/2, 4/2) and Ff compatible with Z € so that (63) holds and ||vf - vf || < 64 < e. 

A.2 Proof of Theorem [3] 

Since Z is SLASF, then it is 7-SLASF for some < 7 < 1 . Then, Theorem[TJshows that Z £ is 7-CEAS 
for I = 1 , . . . , n, irrespective of which and F^ 7-compatible with Z l are taken at each iteration L Since 
Z is controllable, then (A, 1 = A,, = B,) is controllable and (A ( £ , B ( £ ) is controllable for £ = 1 , . . . , n (see, 
for example, Proposition 1 .2 of uM). 

Since (Aj, Bj) is controllable and B ( - g C nx1 , then B ( - ^ and B, has full column rank. By the continuity 
of ©, given e > we can find ei > such that ((36) holds whenever 

||A-A||<£1, IIB--B/H <ei, (64) 
\\Ki- K,\\ <e v (65) 

By ||8), given e-\ > we can find < e 2 < e-\ such that (65)1 holds whenever 

\\Ff - Ff\\ < e 2 , \\0 e - U e \\ < e 2 , \or £ = ^, . . . , n. (66) 

For £ = n, A? and Bf 1 are scalars and Bf ^ because (/A ( ", Bf) is controllable. Consider the following 
condition 

HA n -A"ll <ea. I|e, n -e, n || <e 3 - (67) 

Note that if we select < e 3 < e 2 small enough then {67]| will imply that Bf ^ 0. Taking such small e 3 , 
for each set {(A/ 1 , B ( ") : / G n} satisfying d67j Procedure CEA in Algorithm [TJ does not employ e n d , and 
returns Ff such that |A"' CL | < 1 - e", which is possible for every < e" < 1. Note that if we choose 
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e 3 > small enough, Ff will, in addition, satisfy \A" + BfFf 1 ] < 1 - e n c /2. Therefore, we may take 
F? = F». 

If n = 1,0 is established taking e* = 1, arbitrary e* > 0, and 8 = e 3 , since in this case K, = Ff and by 
the above argument we may take K, = K,. 

We proceed for n > 1 . We next establish the existence of e* and e*. For £ = 1 , . . . , n — 1 , i.e. 
for n r > 1 , consider the sets Tjf as defined in ((27). Note that the set Xjf depends on the matrices Af 
and Bf, for / e n. Since the latter matrices depend on the specific vectors and matrices computed at 
previous iterations of Algorithm [Q then "Ti also depends on such quantities. For £ = 1, . . . , n — 1 and 
k = 1 , . . . , £, consider the expression 



rfr' 1 := inf inf 



inf min d(vf, Im Bf) 



(68) 



Note that the expression between square brackets in {68) coincides with e e f as given in ((29). The 
expressions (1681 can be interpreted as the infimum of the minimum distance between all the possible 7- 
compatible vectors that could be obtained at iteration £ and all the possible Im Bf that can be obtained at 
iteration £, having the data corresponding to iteration k < £ and over all possible outcomes at iterations 
k,k + \,...,£. We make the following claim, to be proved later. 

Claim 4. For I = 1 , . . . , n — 1 and k = \,...,l, expression {68|) 

a) may depend on Af, Bf and 7 but does not depend on Ff or for k < r < £ - 1 . 

b) is positive. 

By Claim HI knowing Af = A,, Bf = B-, and 7, we can define positive 

4 ■= 7, (69) 

e -=, =1 m 5-i^- (70) 

Note that e* < e £ f for £ = 1 , . . . , n - 1 , with e e f as in (29). 

Applying Theorem [2] repeatedly (n - 1 times), it follows that each < e£ < e* and < e e d < e* d 
ensure that given 5 n := £3 > 0, we can find a corresponding < 8g. < 63, for £ = n — 1 , . . . , 1 so that for 
each Af, Bf satisfying 

\\Af - Af\\ < 5 e , \\Bf - Bf\\ < 5 e , (71) 

there exisQ \/\ e S e (e e c /2, e^/2) and Ff compatible with Z e , and a unitary matrix {26) that cause, for 
£ = \,...,n-\, 

\\Ff -Ff || <6 M> ||0f-vf||<fc +1) (72) 
|| U M - U M || < 8 M , \\Af +: - Af + ' || < 5 M , ||Bf +1 - Bf +1 1| < . (73) 

By definition, U-\ = U-[ = I. The latter fact jointly with ((72> -J73t establish (HD. Take <5 = 8-\. Since 
(/f E 5 £ (e^/2, e e d /2) and vf , Ff are compatible with then Lemma|3]shows that {Af + BfFf)vf = Afuf 



1 We add a superscript I to the set 5 to denote that this set is different at each iteration of Algorithm [TJ 
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with \\* | < 1 - e£/2, for £ = 1 n — 1 . Also, is an eigenvalue of /A ( 0L (recall RemarkQ), and so is 

A? + B"F", which satisfies \AP + B*?Pf\ < 1 - eg/2. Therefore, p(/A ( 0L ) < 1 - e c /2 and the proof of partQ) 
is concluded if Claim H] is true. 

To establish ClaimH] consider the sets Xjf defined in (|2S), for I = 1 , . . . , n — 1 . Note that 7^f depends 
on >Af, Bf which, by and (l9l-(fT2t. may depend on A^~\ B l ~\ P^ and LV We next establish that 
7^f does not depend on Ff~ 1 . From 0, CLD and CLD, we have 4/ = U*Af~^ U e + BfPf~ ] U e . Consider 
Pf as in ((43J, and note that PfBf = 0. Thus, it follows that Pf(\\ - Af) = Pf(X\ - U*Af^U e ), which 
does not depend on F^'" 1 . From the equivalent definition of Xjf in (|42 | -(|43 ) , it follows that Zjf does not 
depend on F^~ 1 . Following similar considerations, we can show that Xjf does not depend on Pf for 
k = \,...,t-\. 

Next, consider r)if from {68}, rewritten to explicitly show the dependence on matrices: 
iB = inf mind(v £ , lmL/;Bf^ 1 ), (74) 

We next show that for fixed A l ~ A and B ( £ ~ 1 , ((74) depends only on the orthogonal complement of Ue and 
not on U e itself. Consider fixed Af~^ , B ( f _1 , Ue, and let be such that [v^" 1 lO^] is unitary. Let Ue also 
make [vf~ 1 |Ly unitary. Write B ( £ ~ 1 as 

Bf- ] = vf~ 1 a + L/,/3, = vf " 1 a + Cltf, (75) 
Using ((75) we may write 

d(v e , Im O^Bf- 1 ) = d(^, Im A) = 6(v e , Im 0;0^ ; ) (76) 
Put v e = U*UeV l and, since U*Ue is unitary, we have 

d(v e , Im /3,) = d(v e , Im /?,), and min d(^, Im /3,) = min d(i^, Im /3,). (77) 

/Gn i'6n 

Let = 7^(/4f- 1 , Bf- : , Q*) and note that v l e ff if and only if v e e 7?. Therefore, 

inf mind(v £ , Im /?,) = inf min d(v e , Im /3,) (78) 

Recalling that /3; = U* t B^ and /3, = U^Bf~\ we establish that (74) depends on 0^ only through its 
orthogonal complement Following similar considerations, we can establish that rfcf depends on 

U r only through its orthogonal complement v r ^ , for r = k t. 

To prove part© of Claim consider (ZD again. Since r}if = t l £ as in we already know that 
r\if > for every possible . Next, for fixed /A^~ 1 , B £ ~ 1 , we may write 

rtf = -g(vf" 1 ) = inf -f(0,,vf) = - sup f(0,,vf), (79) 
rfr u = f inf -gtfr\ (80) 
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where, in (75), L/^ must be taken so that [v\ : \Ui\ is unitary, and we have defined 

f(0 t , vf) = - mind(vf, Im i^Bf- 1 ), (81) 



K(U e ) = % e (Af-\Bf~\U e ). (82) 

Note that f is continuous at every (L\ v £ ) for which 0*6^ 1 * for / G n. 
Let W denote the following set 

^ := {U e G c n ~^ n ~^ 1 : [vf -1 is unitary, vf" 1 G Xf~ 1 } (83) 

and regard TZ as the set-valued map TZ:U ^ S-\ that maps each G U to the set 7£(0^) c Si . Making 
these considerations, the function a : 7^~ 1 ->■ R defined above can be regarded as a marginal function. 
Application of Theorem 1.4.16 of JTsfl to g yields that g is upper semicontinuous if f and TZ are upper 
semicontinuous, and TZ has compact values. Since f is continuous, it is upper semicontinuous. TZ takes 
compact values because Xjf is compact. 

Claim 5. The set-valued map TZ : U ~» <Si /s upper semicontinuous. 

Proof. Since 7£ has compact values, i.e. the sets 7^(0^) are compact for every 0^ G W, we need only 
prove that the graph of TZ is closed. For, take a sequence {u k G ^}^ so that \\m k ^, 00 u k = Ue and 
suppose that v k g 7£(u/<) so that \\m k ^, 00 v k = v g <Si. We have to establish that v g TZ(Ue). Recalling 
(|42 ) - i|43) , we have that v k G TZ(u k ) implies 

Pl k (\ iM \ - u* k A{-' u k )v k =0, (84) 

with |A/,/c| < 1 — 7 and 

^ = ['-^f- 1 (K6f- 1 )*u^/- 1 f 1 (u^ 
for all / G n and k > 0. We may rewrite ((84) as 

Pf ik v k \ i>k = Pf k u* k Af-'u k v k . (86) 
Since U/< and v k are convergent, the limit of the right-hand side of ((86) exists and we have 

lim Pf k v k \ i>k = PfU*Af-'U i v= lim Pfv\ k , (87) 

where the last equality above follows because \j tk is bounded. If Pfv + 0, then lim^oo X i:k = A,, with 
| A,| < 1 - 7. If Pfv = 0, then ((87) equals zero. In either case, we have v g TZ(Ue) = f*(Af-' 1 , Bf _1 , 0^). 

□ 

The continuity of f and Claim [5] imply that g is upper semicontinuous. Then, —g is lower semicon- 
tinuous and since T^~ A is compact, then -g will attain a minimum within 7^f _1 . Therefore, > 0. 
Following similar considerations, we can establish that r^ l ' t > for k = 1 , . . . , 1. 

O By Definition |U a SLASF Z with compatible K, means that Af L will generate a solvable Lie algebra 
and satisfy p{A^) < 1 . In addition, partE) ensures that p(Af l ) < 1 - e c /2 irrespective of which e or 5 are 
selected. By LemmaQ] the closed-loop DTSS with subsystem matrices Af L admits a CQLF PartE) then 
follows from robustness of the CQLF, selecting e > small enough. 



?^~ 1 )* 



(85) 
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